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Abstract 

We prove the Gallavotti-Cohen fluctuation theorem for a model of heat 
conduction through a chain of anharmonic osciUators coupled to two Hamil- 
tonian reservoirs at different temperatures. 



1 Introduction 

The Gallavotti-Cohen fluctuation theorem refers to a symmetry in the fluctuations 
of the entropy production in nonequilibrium statistical mechanics. It was first 
discovered in numerical experiments of Evans, Cohen and Morris [|| and then 
discussed in |^ in the context of thermostated systems. As a mathematical theorem 
it was proved for Anosov dynamical systems . Soon thereafter the fluctuation 

theorem was discussed in the context of stochastic dynamical systems first by 
Kurchan lIT^ and then, more systematically by Lebowitz and Spohn, and Maes 
[ p2| [l8| . In particular, Maes discovered a general formulation of the fluctuation 
theorem in the context of space-time Gibbs measures which covers both Markovian 
stochastic dynamics and chaotic deterministic dynamics (via a Markov partition). 
As a mathematical theorem the fluctuation theorem is proven for quite general 
stochastic models with finite state space, such as lattices gases in a finite box. 
Relations for the free energy related to the fluctuation theorem have been also 
discussed in jl5[ ||. 

Among the consequences of the fluctuation theorem is the non-negativity of 
entropy production although the proof of its positivity is more difficult and is so far 
proved only in particular examples ^ . We also note that in the related context 
of open systems, classical and quantum, the production of entropy is discussed at 
a general level in |27[ |l^, |24j . Again the non- negativity of entropy production is 
relatively easy to establish, while the strict positivity has been established only in 
particular models 0, p^ . 

In this paper we consider an open system consisting of a finite (but of arbitrary 
size) chain of anharmonic oscillators coupled at its ends only to reservoirs of free 
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phonons at positive and different temperatures [|[ 0, ||, ^ ^ . In particular our 
model is completely Hamiltonian and its phase space is not compact. 

In order to establish the fluctuation theorem, two ingredients are needed: 
one needs to prove a large deviation theorem for the ergodic average of the entropy 
production and establish a symmetry of the large deviation functional. The second 
part is usually relatively straightforward to establish, at a formal level, since it 
follows from a symmetry of the generator of the dynamics. This formal derivation 
for models related to ours can be found in |^ and fl^ . 

The first part, proving the existence of the large deviation functional, involves 
technical difficulties, in particular if the phase space of the model is not compact. 
In this case large deviation theorems are established provided the system satis- 
fies very strong ergodic properties (such as hypercontractivity) see e.g. |§, ||, p9| . 
In addition the entropy production is in general an unbounded observable while 
standard results of large deviations apply only to bounded observables. 

In this paper we show how to treat these difficulties in the model at hand. 
The techniques we use are based on the construction of Liapunov functions for 
certain Feynman-Kac semigroups and Perron-Frobenius-like theorem in Banach 
spaces. We heavily rely on the strong ergodic properties of our model established 
in 1^, 1^, ||] and especially in |^ . 

The Hamiltonian of the model, as in Q , has the form 

H = Hb+Hs + Hi. (1) 

The two reservoirs of free phonons are described by wave equations in R"* with 
Hamiltonian 

Hb = H{ipl,ttl) + H{ipR,TrR) , 
Hi^^n) = iy'da;(|V^(x)|2 + Kx)p), 

where L and R stand for the "left" and "right" reservoirs, respectively. The Hamil- 
tonian describing the chain of length n is given by 

n 2 
n n — 1 



i=l 



where {pi,qi) E R'' x R'' arc the coordinates and momenta of the particle of 
the chain. The phase space of the chain is R^''". The interaction between the chain 
and the reservoirs occurs at the boundaries only and is of dipole-type 

Hi = qi ■ J dxV(pL{x)pL{x) + Qn ■ J dxVipR{x)pR{x) , 
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where pL and pR are coupling functions ("charge densities"). 

Our assumptions on the anharmonic lattice described by Hs{p,q) are the 
following: 

• HI Growth at infinity: The potentials C/^^) (x) and U'-'^^ (x) are C°° and 
grow at infinity like WxW''^ and HxH*^^: There exist constants Ai, Bi, and Cj, 
i = 1,2 such that 

lim \''''U^'^Xx) = A,\\xf^ , 

A — >-oo 

lim A-'=-+iV[/W(Aa;) = Aiki\\xf'-'^x, 

A— j-oo 

||a2c/W(x)|| < {Bi + CiVix))'-^i . 
Moreover we will assume that 

fc2 > /Cl > 2, 

so that, for large ||x|| the interaction potential C/^^^ is "stiffer" than the 
one-body potential U^^\ 

• H2 Non-degeneracy: The coupling potential between nearest neighbors 
J/^^^ is non-dcgcncratc: For x G R'' and m = 1,2, let A(™)(x) : R'' 
R''™ denote the linear maps given by 

We assume that for each x G R"^ there exists mo such that 

Rank(A(i)(x),---A('"«)(x)) =d. 



• H3 Rationality of the coupling: Let pi denote the Fourier transform of 
Pi. We assume that 



where Qi, i e {L,R} are polynomials with real coefficients and no roots on 
the real axis. 

We introduce now the temperatures of the reservoirs by choosing initial con- 
ditions for the reservoirs. The Hamiltonian of a reservoir is quadratic in = (cp, tt), 
H = 5')/2, and therefore the Gibbs measure at temperature T, dpT{'^) is the 
Gaussian measure with covariance T {■,■). To construct nonequilibrium steady 
states we assume that 
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• The initial conditions '^l — {4>l,''^l) and "^r = {(/)r,ttr) of the reservoirs 
are distributed according the gaussian Gibbs measures dfiTL and d/XTj, re- 
spectively. 

In order to define the heat flow through the bulk of the crystal we consider 
the energy of the z*'' oscillator which we take to be 

H, = § + U'^^\qi) + \ (c/(2)(q,_i - q^) + U^^\q, ~ g.+i)) . (2) 

Differentiating Hi with respect to time, one finds that 

dH, 



where 



dt 

$. = ^^^^^V[/(^)(g. (3) 

is the heat flow from the i*'* to the (i + 1)*'* particle. We define a corresponding 
entropy production by 

where Tr and Tl are the temperatures of the reservoirs. 

There are other possible definitions of heat flows and corresponding entropy 
production that one might want to consider. One might, for example, consider the 
flows $1,, ^R at the boundary of the chains, and define ai, = —^l/Tl — ^r/Tr, 
or one might take other quantities as local energies. But using conservation laws 
it is easy to see that all these heat flows have the same average in the steady 
state. Moreover we will show that all the entropy productions have the same large 
deviations functionals: the exponential part of their fluctuations are identical. 

We denote {p{t),q{t)) = (p(i,po, 9o, 9(i,Po, 9o, ^'l, as the Ha- 

miltonian flow generated by the Hamiltonian (p, and consider the ergodic average 

1 ft 

(Ti' = 7 / cr^{p{s),q{s))ds . 
^ Jo 

The quantity ai(j){s),q{s)) depends on both the initial conditions of the chain 
and of the reservoirs which, by assumption, are distributed according to thermal 
equilibrium. By the ergodic theorem proven in there exists a measure dv on 
R^'*" such that 

lim ct/ = I Uidv . 



for all (poi'Zo) and d^,T^ and d^TR almost surely. Moreover Jaidv = (ct)^ is 
independent of i and as shown in j7| 

(cr)^ > and (o-)^ = if and only if Tl=Tr. 
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Given a set A C R, we say that the fluctuations of Ui in A satisfy the large 
deviation principle with large deviation functional I{w) provided 

inf I{w) < liminf-ilogP{CT^* e A} < 

■u;eInt(A) t^oo t 

limsup-ilogP{CTj* e yl} < inf liw). 

t^oo t tuGClos(A) 

The study of large deviations for ai is based on the moment generating functionals 
ei{a) given by 

e.(a) = lim -^log / d^T.^MT^e^/o ^'(^(^^'^^^^ 

t^cc t J 

The main technical result of this paper is 
Theorem 1.1 // 

_ I ^min ^rr 

a e — , 1 



max mill max min 

e{a) = ei{a) is finite and independent of i and the initial conditions {pQ,qo). 
Moreover e{a) satisfies the relation 

e{a) = e(l — a) . 

As an application of the Gartner-Ellis Theorem, see [Q, Theorem 2.3.6, we 
obtain the Gallavotti-Cohen fluctuation theorem. 

Theorem 1.2 There is Q, Ticighbovhood O of the intcTval ; 

1 such that 

for A d O the fluctuations of ai in A satisfy the large deviation principle with a 
large deviation functional I{w) obeying 

I{w) — I{~w) = —w , 

i.e., the odd part of I is linear with slope —1/2. 



Theorem 1.2 provides information on the ratio of the probabilities of observing 



the entropy production to be w and —w: roughly speaking we have 



P{a,* e (w-e,u; + e)} _ 



P{ct/ e {-w - e, -w + e)} 

2 Fluctuations of the entropy production 
2.1 Exponential mixing and compactness 

As shown in ^, ^ , under condition H3 the dynamics of the complete system can 
be reduced to a Markov process on the extended phase space consisting of the 
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phase space of the chain and of a finite number of auxihary variables which 

we denote as r. In the simplest case which we consider here, (corresponding to 
p{k) ~ (fc^ + 7^)~^)7 f = {i"!, ''ri) <= R^'^ and the resulting equations of motion take 
the form 

q = P, 

p = „v,F-A^r, 

dr = {-jr + Ap)dt+{2-/Ty/^duj. (4) 

Herep = {pi, ■ ■ ■ ,pn) and q — (gi, • • • , g„) denote the momenta and positions of the 
particle, r = (ri, Vn) are the auxiliary variables and w is a standard 2(i-dimcnsional 
Wiener process. The linear map A : R''" R^'' is given by A(pi,...,p„) = 
{Xpi,Xpn) and T : R^'^ ^ R^'' by T{x,y) = {Tix,T„y). Here Ti = Tl and 
T„ = Tr are the temperatures of the reservoirs attached to the first and n*^ 
particles respectively, 7 is the constant appearing in p and A is a coupling constant 
equal to HpHls. 

The solution of Eq. (|), x{t) = {p{t), q{t),r{t)) with x e X ^ j^2d(«+i) ^ 
Markov process. We denote T* as the corresponding semigroup 

TV(x) = E,[/(xW], 

with generator 

L = 7 {WrTWr ~ rWr) + (ApV, - rAVp) + (pV, - (V,F(g))Vp) , (5) 

and we denote Pt{x, dy) as the transition probability of the Markov process x{t). In 
[ p6| we proved that the Markov process x{t) has smooth transition probabilities, in 
particular it is strong Feller^ and that it is (small-time) irreducible: For any t > 0, 
any x & X and any open set A C X we have Pt{x, A) > 0. 

There is a natural energy function associated to Eq.(^, given by 

G{p,q,r) = — + H{p,q), 

which we employ throughout our discussion. In [p6| we have constructed a Lia- 
punov function for x{t) from G: Let t > and < 9 < max(ri, T„)^^. Then 
there exists Eq such that for all E > Eq there exist functions k = k{E) < 1 and 
b = b{E) < 00 such that 

T'e'^'ix) < K{E)e'''ix) + b{E)l[G<E}ix) . (6) 

Moreover k(E) can be made arbitrarily small by choosing E sufficiently large, in 
fact there exist positive constants ci = ci{9,t) and C2 = 02(0, t) such that 

k{E) < ce-^^'^'^'r (7) 
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By results of it is also shown in that the convergence to the unique 
stationary state, denoted by /i, occurs exponentially fast: Let Tioofi denote the 
Banach space {/; ||/||oo,e = sup^ \f{x)\e^'^'^''^^ < oo}. Then there exist constants 
r > 1 and R < oo 



\T* fix) - I fd^x\ < Rr- 



l/ll 



(8) 



which means that T*, acting on 7ioo,e has a spectral gap. The methods of are 
probabilistic and rely on a nice probabilistic construction called splitting as well 
as coupling arguments and renewal theory. 

Under the condition given here, by taking advantage of the fact that the 
constant n in the Liapunov bound can be made arbitrarily small (this is not 
assumed in we can prove stronger ergodic properties and also give a direct 

analytical proof of Eq. (||) . 

Besides the Banach space TLoo,e defined above we also consider the Banach 



space Vy^ 



{/, |/|e-«G G Co(X)} with norm 



( Cq{X) denotes the set 



of continuous functions which vanish at infinity). Furthermore for 1 < p < oo we 
consider the family of Banach spaces Tipfi = LP{X, e^^^'^'^^^dx) and denote || • \\pfi 
the corresponding norms. 

Theorem 2.1 IfO< OT^ < 1, the semigroup T* extends to a strongly continuous 
quasi-bounded semigroup on Hp^g, for 1 < p < oo and on Ti.'^ g. For any t > 0, 
is compact on Tipfi, for \ < p < oo and on H^^ g. 

As an immediate consequence of the spectral properties of positive semi- 
groups and the irreducibility of x{t) we have 

Corollary 2.2 The Markov process x{t) has a unique invariant measure dp, and 
Eq. d; holds. 

Proof: : Since T* is a Markovian, compact, and irreducible semigroup the eigenvalue 
1 is simple with the constant as the eigenfunction. This shows that the Markov 
process x(t) has a unique invariant measure. Moreover by the cyclicity properties 
of the spectrum of a positive semigroup and by the compactness of T*, there 
are no other eigenvalues of modulus 1. Eq. M) follows immediately. | 



Proof of Theorem |g. J| In Ig^li Lemma 3.6, we sh owed that for some constant C 
rpt^ec < ^ct^ea provided OT, < 1 (see also Lemma |2j| below). Therefore for / C°° 
with compact support we have, using Ito's and Girsanov's formulas 



-^GT*e««/(x) 



J{G(x(t))~G(x)) 



^'feTi{T)+'ff{e^T-9)f 



fim) 
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where x is the process with generator 

Le = L + 2j0rTWr ■ 

A computation shows that Ljl = 7Tr(l — 26T). Standard arguments show then 
that the semigroup associated with the process x extends to a quasi-bounded 
and strongly continuous semigroup on L^[dx), 1 < p < oo and on Cq{X). Using 
the assumption that 9Ti < 1 and Feynman-Kac formula we see that e^^'^T^e^'^ 
extends too to a quasi-bounded and strongly continuous semigroup on LP{dx), 
1 < p < oo and on Co{X). This implies immediately that T* extends to a strongly 
continuous semigroup on Hp.e, I < p < oo and g. The computation above also 
shows that T* extends to a quasi-bounded semigroup on TLoo.e- 

We first prove the compactness of T* for Hoo.e- If / G T^oofi then \ f{x)\ < 
||/||oo,ee^'^(^) and by (|) and §) we obtain 

\lG>ET'm\ < e««(^) sup ^!JM 

{y:G(y}>E} e'^'-^V) 

rpt (0G(y)) 

< e«G^^^II/l!oo,. sup I 

{y:G{y)>E} e«G(y) 

< K{E)e''"^^^\\f\Ue. (9) 

From the bounds (^) and we conclude that the operator 1{g>e}T'* converges 
uniformly to in 7ioo,e as ^ oo. The semigroup T* has a C°° kernel since it is 
generated by a hypoelliptic operator see Proposition 4.1, so, by the Arzela- 
Ascoli theorem 'i-{G<E}T^^^^{G<E} is compact, for any E. Therefore we obtain 

T* = hm 1[G<E}T'/H[G<E}T'^^ , 

E — >oc 

where the limit is in the norm sense from above, i.e., T* is the uniform limit 
of compact operators, hence is compact. 

The compactness of T* for g follows from the same argument. In fact by 
Eq.®, for any t > 0, T'Hoofi C ^S^^g. 

To prove the compactness of T* on Tipg , 1 < p < oo , we note that 

|r*/(x)| = \EM(.^m\ 

= |E,[ef«(^(*»e-f«(-(*))/(x(t))]| 



Thus using the bound and the fact that T* is quasi-bounded on Hi^ we obtain 
\\lG>ET'mr> < I E.[e««(-W)]f E.[e-f«(-(*»F(a:(t))]e^^'^«(-)da; 

J {x:G(x)>E] 
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/E fe'"='(^(*))l\ ' per. 

{:e:G(x)>B} V e^'^^'"^ / 

< «(ii;)fe^*||e-^«/f||i,e 

As in the case p — oo, we conclude from the bound that the operator 
^g>eT^ converges uniformly to in TLp^ as E ^ oo. Using that the kernel of 
^{G<E}T*^{G<E} is bounded, we conclude that T* is compact on Tip^ for 1 < p < 
oo. I 



2.2 Heat flow and generating functionals 

In order to define the heat flows we note that we have 

j^T'H = LT'H = T\~rAp) = T*(-Aripi - Ar„p„) . 

Hence we identify <i>o = — Aripi as the observable describing the heat flow from 
the left reservoir into the chain and $„ = Ar„p„ as the heat flow from the chain 
into the right reservoir. As in the introduction we define the energy Hi of the i*"^ 
oscillators by Eq.(^, for i < 2 < rt — 1, and 

H, ^ | + c/(i)(,0 + lc/(2)(,,_,,), 

With the heat flows i = 1, • • • , n, defined as in Eq. (^) we have 

LH, = - <I>j , i 1, • • • 
and we define the entropy productions (7^, i = 0, • • • , n by 

a.= (^-^)$. . = 0,...,n. 

We now provide several identities involving the generator of the dynamics and 
the entropy production, which will play a crucial role in our subsequent analysis. 

Lemma 2.3 Let the function Ri, i = 0, ■ ■ ■ ,n be given by 

^^ = ^(?+E^'(P''?))+^( E H,ip,q) + Yj ■ (10) 
\ k—1 / z+1 / 

Then we have 

ai^rT-^r-TT{I)+LRi. (11) 
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Proof: This is a straightforward computation. | 



Remark 2.4 This shows that, up to a derivative, aU the entropy productions are 
equal to the quantity rT^^r — Ti^I which is independent of i and involves only 
the r-variables. 

Let be the formal adjoint of the operator L given by Eq. (H) 

= 7 {VrTVr + rVr) " (ApV. - rAVp) - (pV, - (VgF(g))Vp) , (12) 

and let J be the time reversal operator which changes the sign of the momenta of 
all particles, Jf{p,q,r) = f{-p,q,r). 

The following identities can be regarded as operator identities on C°° func- 
tions. That the left and right side of Eq. (|l^) actually generate semigroups for some 
non trivial domain of a is a non trivial result which wc will discuss in Section 2.3. 

Lemma 2.5 We have the operator identities 

e^'JL'^Je-"' ^ L-a,, (13) 

and also for any constant a 

e-""J{L^ - aa^)Je^^ = L - (1 ~ a)a, . (14) 

Proof: We write the generator L as L = Lq + Li with 

Lo = 7(VrTV,r - rVr) (15) 
ii = {ApVr - rAVp) + [pVq - {VqV{q))Vp) . (16) 

Since Li is a first order differential operator we have 

e"^'Lie^' = Li+ LiRi = Li + a.i . 

Using that WrRi = T~^r we obtain 



e 



This gives 



= 7V^r(V^+T-V) = L'o . 



which is Eq. (p^). Since JuiJ — —ai, Eq. ( p4| ) follows immediately from Eq. (13). 
Remark 2.6 In the equilibrium situation, i.e., for Ti = T„ = T, Eq. (|lj) is 
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which is simply detailed balance. Eq. (Q) can be interpreted in path space in the 
following manner jl8|: Let 11 denote the time-reversal in path space on the time 
interval [0,t]: U{p{s), q{s),r{s)) = {—p{t — s),q{t — s),r{t — s)) and let dP denote 
the measure on C([0,t],X) induced by x{t). Then Eq. ( p4[ ) implies that 

dPoU _ R,{x{t))-Ri{x{Qy)-J*criixis))ds 



dP 



This formula exhibits the fact that the lack of microscopic reversibility is intimately 
related to the entropy production. 

We now turn to the study of the large deviations. As shown in the Markov 
process x(t) is ergodic. In order to study the large deviations of t^^ ai{x{s))ds 
we consider the moment generating functionals 



rut, a) = E, 



Formally the Feynman-Kac formula gives smce (Ti 

is not bounded, nor even relatively bounded by L, it is not obvious that Tl.(t,a) 
exists for a Q. Our goal is to prove that r^(i, a) exists and that the limit 

e(a)= hm -ilogr^(i,a) (17) 

exists and is finite in a neighborhood of the interval [0, 1], and is independent of i 
and of the initial condition x. 

The technical difficulty in proving the existence of the limit ( p7| ) lies in the 
fact that the functions ai are unbounded. Standard large deviation theorems for 
Markov processes (see e.g. ^, ^, ^) are proven usually under strong ergodic 
properties for bounded functions and are not directly applicable. Large deviations 
for unbounded functions are considered in 0] for discrete time countable state 
space Markov chains under conditions which amount in our case to cr = o(G). In 
our case this is clearly not satisfied since, in general a is not bounded by G. 

But the ai are very special observables, in particular they are intimately 
linked with the dynamics as shown by the identities Eqs.(|3|) and (|lj). The next 
lemma displays another identity which will be important in our analysis. 

Lemma 2.7 We have the identity 

L-aa, = e"^'L„e~"-"' , (18) 

where 

la = - ((a - a^)jrT-^r - aTr(7/)) (19) 

and 

La = L + 2ajrVr. (20) 
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Proof: As in Lemma 2.5 we write the generator L as L ~ Lq + Li, see Eqs.(|T^) 
and (|l|). Since Li is a first order differential operator we fiave 

e-"^'Li_e"^' = Li+a{LiRi) = L^+aUi. (21) 

Using that VrRi = T^^r is independent of i we find that 

— Lq + a7(rVr + V^r) + (a^ — a)^rT^^r 

= Lo + 2Q;7rVr + (a^ - a)7rr" V + Q;Tr7/ . (22) 

Combining Eqs. and (p2[) gives the desired resuh. | 



Remark 2.8 The identity (y_8|) shows that all operators L — aai are conjugate 
to the same operator La- This will be the key element to prove that e(a) is 
independent of i. Furthermore it can be seen from Eqs. (^9|) and ( pO[ ) that La has 
the form of L plus a perturbation which is a quadratic form in r and Vr-. Such a 
perturbation is indeed nicer than aai. Also it should be noted that La has very 
much the same form as the operator L: they differ only by the coefficient in front 
of the term rV^. This fact will allow us to use several results on L obtained in 

2.3 Liapunov Function for Feynman-Kac Semigroups 

At this point we begin the study of La as the generator of a semigroup. 

Proposition 2.9 If 6 and a satisfy the condition 

-a<6T,<l-a, (23) 

then there exists a constant C — C{a,6) such that e*"^°e^*^(x) < 6*^*6^*^(2;). 

Proof: We note first that La, defined in Eq. (p^), for all a S R, is the generator 
of a Markov process which we denote as x{t). Indeed we have that 

LaG{x) = Tr(7T) - (1 + 2ay < Ci + C2Gix) 

Since G grows at infinity, G is a Liapunov function for x{t) and a standard argu- 



ment 1 16 shows that the Markov process x(t) is non-explosive. Furthermore we 



have the bound 



Laexp9G{x) = 

= exp 0G{x)-f [Tr(6'r + al) + r{e^T - (1 - 2a)0 ~ a{l - a)T-'^)r] 

< CexpOGix) , (24) 
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provided a and Ti, i — l,n satisfy the inequality 



e^T, - (1 - 2a)0 - a(l - a)Tr^ < , 



or 



<9T, <l-^ 



We denote cr^ as the exit time from the set {G{x) < R}, i.e., aji = inf{t > 
0,G{x{t)) > R}. If the initial condition x satisfies G{x) = E < R, we denote by 
Xji{t) the process which is stopped when it exits {G{x) < R}, i.e., Xji{t) = x{t) 
for t < aji and iij(t) — x{aif ) for t > gr- Finally we set <Tii{t) = minja^j, t}. 

By Eq. (^), the function W{t,x) = g-ctg^Gi^) satisfies the inequality {dt + 
La)W{t, x) < and applying Ito's formula with stopping time to the function 
W{t, x) we obtain 



E 

and thus 



(("-a')7fT-if-aTr(7/)) ds 0G{x{an{t))) 



Since the Markov process x{t) is non-explosive G{xii{t)) 
as i? — > oo, so by the Fatou lemma we have 



G{x{t)) almost surely 



This concludes the proof of Lemma E^. | 



The next three theorems are all consequences of the fact that is the 
generator of a Markov process which is similar to the process generated by L: 
Indeed L and Lq differ only by the coefficient in front of the rV^ term. Therefore 
repeating the proofs of |^ we obtain 

Theorem 2.10 The semigroup e*^° has a smooth kernel qait,x,y) which belongs 
to C°°((0,oo) X X X X). 

Proof: The operator La satisfies the same Hormander-type condition that L proven 
in 111, Proposition 4.1. The result follows then from ||l| or H]. | 



Theorem 2.11 The semigroup e*^° is positivity improving for all t > 0. 

Proof: The semigroup e*^° is shown to be irreducible exactly as e*^, see [0, 
using explicit computation and the Support Theorem of ||2^ . The statement follows 
then from the Feynman-Kac formula. | 

As is apparent from the form of La we will need estimates on the observable 
in the sequel. Such estimates were also crucial in for the construction of a 
Liapunov function. 
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Theorem 2.12 Let < a < l/k2 and let ts = £;i/fe2-i/2^ y^^g^g g^-j^^^ ^ ^g^ 

Six, tE) C {/ e C([0, tE],X) ■ /(O) = a:, G{x) = £;} , 
a77.rf constants Eq < oo and A, B , C > such that for E > Eq 

P{xeS{x,E,tE)} > 1 _ ^6-^^'°^'''"'"% 

and 

l-tE 

r^{s)ds > CE^/'''-^/^ , if ie S{x,E,tE). (25) 



Proof: The proof is exactly as in ||26|] . One first sets Ti — Tn — in the equations of 
motion and then, by a scaling argument, Theorem 3.3 of one shows that the 
deterministic trajectory satisfies the estimate (p5|). Then one shows, see Proposi- 
tion 3.7 and Corollary 3.8 of that the overwhelming majority of the random 
trajectories follows very closely the deterministic ones. We refer the reader to |2^ 
for further details. I 



Remark 2.13 For large energy E, paths satisfying the bound ( |25| ) have a very 
high probability. From Eq. ( |25|) we obtain that, on a time interval of order 1, 

f2(s) > CE'^/''' , 

for an overwhelming majority of the paths. 







Theorem 2.14 Let t > be fixed and suppose that a and 9 satisfy the condition 
Eq. (p^. There exist a constant Eq and functions k(E) and b{E) such that for 
E> Eo 

e'^-e'^{x) < K{E)e'^^^^ + b{E)l^a<E}ix) . (26) 
Moreover there exist constants c\ and C2 such that 



k{E) < cie 



Proof: By Propositi on |2.9| the function e^^''e^'^ [x) is bounded on any compact set. 
Therefore to show (p6[) it suffices to show that 



sup 

{x:G{x)>E'} 



' J*{a{l~ahfT 'r-aTr(7/))ds fl(G(£(t))-G(£)) 



< k{E). 



Using Ito's formula we have 



G{x{t)) - G(x) = / 7(Tr(T) - f^) ds + / ./2^rduj{s) , 
Jo Jo 
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and thus we obtain 
E 



e~ lo f^("(l~")7T-l-7e(l-2a))f ds ^ 



„t7Tr(eT+a/; 



E 



xe 



(27) 



Using the Holder's inequality we find that the expectation on the r.h.s of Eq. (|2 
can be estimated by 



E, 



1/9 



xE, 



E.T 



-4^ jl[y/^f?ds^p /Je(27T)i/^fdc.) 
-q-y J* 7=(a(l-Q)T-i-e(l-2Q)+p0^T)f ds 



1/p 



1/9 



(28) 



where we have used that the second factor is the expectation of a martingale with 
expectation 1. 

If 9 and a satisfy the condition (p3|), then, by choosing p sufficiently close 
to 1, the quadratic form in the right side of Eq. ( p8| ) is negative definite. Using 
Theorem ^.12 as in Theorem 3.11 of we obtain 



sup E^ 



- Jg (a(l-a)rT ^ f-aTv{-f I)) ds 0(G(x(t))-G(x)) 



^7Tr(eT+Q/) -CB^/'"2^Tr(Q(l-Q)T"^-(l-2a)6l+pe^T) 



< cie 



and this concludes the proof of Theorem 2.14 



As in Theorem 2.1 we obtain 



Theorem 2.15 If a and 6 satisfy the condition Eq. (pi^j, then e*^° extends to a 
strongly continuous quasi-hounded semigroup on Ti-pfi for 1 < p < oo and on Tl'^ g. 

Moreover e*^" is compact on Tip^g, \ < p < oo and on Ti.^ g. 



Proof: The proof is a repetition of the proof of Theorem 2.1 and is left to the 
reader. | 



As a consequence of Theorem 2.15 and of the theory of semigroup of positive 
operators |ll| we obtain 
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Theorem 2.16 // 



a e 



T ■ T ■ 

^ mm -, I min 



then 



Ti rf~\ ' rr\ rri 

max -'min -'max rr 



e(a) = lim ^^\ogTl{t,a) 



exists, is finite and independent both of i and x. 



Proof: By Theorem 2.15, e*-'"° generates a strongly continuous semigroup on THP^ g 
if 

- a <eT, a. (29) 

If a < 0, this impUes that \a\ < BTa^in < 6'Ti„ax < 1 + |a| and so the set of 6 we 
can choose is non-empty provided 



a > 



If < a < 1, we can always find 6 such that ( |29| ) is satisfied. Finally if a > 1 then 
^ implies that that 

a < 1 



J- max mm 

By the definition of R„ Eq. ® , e'°'"^ £ H^^ g since -a -t- OT, < 0. Using now 
Lemma 2.7, we see that r^(t, a) exists and is given by 



From Theorem 2.11 the semigroup e*^° is an irreducible semigroup of compact 
operators on the Banach space Ti.^ g. From the cyclicity properties of the spectrum 
of irreducible operators and from the compactness it follows (see Chapter 
C-III) that there is exactly one eigenvalue e"**^^"' with maximal modulus and 
this eigenvalue is real and simple. The corresponding eigenfunction fa is strictly 
positive and we denote as Pa the one-dimensional projection on the eigenspace 
spanned by fa- In particular if 5 > 0, then Pag{x) > 0. 

From compactness it follows that the complementary projection (1 — Pa) 
satisfies the bound 



'{l~Pa)fix)\ < Ce-*''(")||/||o,,,e''«(-). (30) 
for some constants C > and d{a) > e{a) and for all t > 0. 



16 



From Lemma 2.7 and Eq. (pO[) we obtain, for all x £ X, that 



lim -]\ogTl{t,a) 

= lim -iloge*(^""'"'h(a;) 



lim — loge"^'e*^°e-"-^'(x) (31) 

i — >oo t 



lim 

t — >oo 



aRi{x) ) + e{a) 



+ lim -- log (p„e-"-"'(x) + e*"(")e*^°(l - P„)e-"^'(x) 

t^oo t \ 

e(a) . 



This concludes the proof of Theorem 2.16 



Using now the identity (|l|) we can prove the symmetry of e{oi). Theorem |l . l| 
is then an immediate consequence of the following result. 



Theorem 2.17 // 



a e 



(32) 



then 



e(a) = e(l — a) . 



Proof: If a is in the interval ( |3^ and —a < 9Ti < 1 — a then e*^" is a strongly 
continuous compact semigroup on TiJ^ g. By Lemma 2.7 



-ct-Ri 



is also a strongly continuous compact semigroup on the Banach space TiJ^ e a ~ 
{/; |/|e-««+"«' e Co(a;)} with the norm ||/||oo,e,a = sup |/|e«G+"«'. 

The dual semigroup (e*(^~"'^»))* is a c ompact semigroup on the Banach space 
(of measures) {H°^^g,J*. By Theorem |t| (e*^^-"'"'))* maps {H^^^g^J* into mea- 
sures with smooth densities and on densities (e*(^^"'^i))* acts as 



(e*(^-""'))*(p(x)da;) = (e 



' p{x))dx . 



By Lemma 2.5 we have 

g-i?.gt(L-(i-a).Oi(a.) - Je*(^"-"^')je-^(x). 



(33) 



Since —a < OTi < 1 — a, e ■'^^ is a density of a measure in {H^^o e a)* ■ Since 
(g*(i-acri)^* -g compact and irreducible with spectral radius e{a) we obtain using 
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Eq. (11) 



e a 



lim --\og\\J{e*^^'^-"''^^Je-^^)dx\\ 

t — *oo t 




sup 



e(l-a). 



In the last equality we have used Theorem 2.16 
measure. This concludes the proof of Theorem 



and the fact that fe ^' is a finite 



2.17 



We finally obtain the Gallavotti-Cohen fluctuation theorem 

Theorem 2.18 There is a neighborhood O of the interval [— (c)^ , (cr)^] such that 
for A G O the fluctuations of at in A satisfy the large deviation principle with a 
large deviation functional I{w) obeying 

I{w) — I(~w) = —w , 

i.e., the odd part of I is linear with slope —1/2. 

Proof: First we note that e{a) is a real analytic function since it is identified with 
an eigenvalue of a compact operator. A simple computation gives that 



a=0 



The function e(a) is analytic and convex. By the result of it is not identically 
zero, and so the symmetry the symmetry e(a) = e(l — a) implies that the set of 
the values of ^e(a) is a neighborhood of [— (cr)^ , (cr)iy]. 

The large deviation principle is a direct application of the Gartner-Ellis the- 
orem, Theorem 2.3.6. The large deviation functional is given by the Legendre 
transform of e(a) and so we have 

I(w) — sup{e(Q!) — aw} — sup{e(l — a) — aw} 

a a 

= sup{e(/3) — (1 — P)w} — I{—w) — w . 

/9 



References 

[1] Balaji, S. and Meyn, S. P.: Multiplicative ergodicity and large deviations for 
an irreducible Markov chain. Stoch. Proc. Appl. 90 123-144 (2000). 



18 



[2] Crooks, G.E.: Path-ensemble averages in systems driven far from equilibrium. 
Phys. Rev. E 61, 2361-2366 (2000) 

[3] Deuschel, J.-D. and Stroock, D.W.: Large deviations. Pure and Applied 
Mathematics 137. Boston: Academic Press, 1989 

[4] Dcmbo, A. and Zeitouni, O.: Large deviations techniques and applications. 
Applications of Mathematics 38. New- York: Springer- Verlag 1998 

[5] Eckmann, J. -P. and Haircr, M.: Non-cquilibrium statistical mechanics of 
strongly anharmonic chains of oscillators. Commun. Math. Phys. 212, 105- 
164 (2000) 

[6] Eckmann, J. -P., Pillet C.-A., and Rey-Bellet, L.: Non-equilibrium statisti- 
cal mechanics of anharmonic chains coupled to two heat baths at different 
temperatures. Commun. Math. Phys. 201, 657 697 (1999) 

[7] Eckmann, J.-P., Pillet, C.-A., and Rey-Bellet, L.: Entropy production in non- 
linear, thermally driven Hamiltonian systems. J. Stat. Phys. 95, 305-331 
(1999) 

[8] Evans, D.J., Cohen, E.G.D., and Morriss, G.P.: Probability of second law 
violation in shearing steady flows. Phys. Rev. Lett. 71, 2401-2404 (1993) 

[9] Gallavotti, G. and Cohen E.G.D.: Dynamical ensembles in stationary states. 
J. Stat. Phys. 80, 931-970 (1995) 

[10] Gentile, G.: Large deviation rule for Anosov flows. Forum Math. 10 89-118 

(1998) 

[11] Grcincr, G: Spectral theory of positive semigroups on Banach lattices. In 
One-param,eter semigroups oj positive operators Lecture Notes in Mathematics 
1184, Ed. R. Nagel, Berlin: Springer, 1986, pp 292-332 

[12] Hormandcr, L.: The Analysis of linear partial differential operators. Vol III, 
Berlin: Springer, 1985 

[13] Jaksic V. and Pillet C.-A.: On entropy production in quantum statistical 
mechanics. Commun. Math. Phys. 217, 285-293 (2001) 

[14] Jaksic V.and Pillet C-A.: Non-equihbrium steady states of finite quantum 
systems coupled to thermal reservoirs. Preprint (2001) 

[15] Jarzynski, C: Hamiltonian derivation of a detailed fluctuation theorem. J. 
Statist. Phys. 98, 77-102 (2000) 

[16] Has'minskii, R.Z.: Stochastic stability of differential equations. Alphen aan 
den Rijn — Germantown: Sijthoff and Noordhofi^, 1980 



19 



[17] Kurchan, J: Fluctuation theorem for stochastic dynamics. J. Phys.A 31, 
3719-3729 (1998) 

[18] Maes, C. The fluctuation theorem as a Gibbs property. J. Stat. Phys. 95 
367-392 (1999) 

[19] Maes, C. Statistical mechanics of entropy production: Gibbsian hypothesis 
and local fluctuations. Preprint (2001) 

[20] Maes, C, Redig, F., and Verschuere, M. No current without heat Preprint 

(2000) 

[21] Meyn, S.P. and Tweedie, R.L.: Markov Chains and Stochastic Stability. Com- 
munication and Control Engineering Series, London: Springer- Verlag London, 
1993 

[22] Lebowitz, J.L. and Spolin, H.: A Gallavotti-Cohen-typc symmetry in the 
large deviation functional for stochastic dynamics. J. Stat. Phys. 95, 333-365 
(1999) 

[23] Norriss, J.: Simplified Malliavin Calculus. In Seminaire de probabilites XX, 
Lectures Note in Math. 1204, Berlin: Springer, 1986, pp. 101-130 

[24] Pillet, C.-A.: Entropy production in classical and quantum systems. Markov 
Proc. Relat. Fields 7, 145-157, (2001). 

[25] Rey-Bellet, L. and Thomas, L.E.: Asymptotic behavior of thermal non- 
equihbrium steady states for a driven chain of anharmonic oscillators. Com- 
mun. Math. Phys. 215, 1 24 (2000) 

[26] Rcy-Bcllct, L. and Thomas, L.E.: Exponential convergence to non-equilibrium 
stationary states in classical statistical mechanics To appear in Commun. 
Math. Phys. 

[27] Ruelle, D.: Entropy production in quantum spin systems. Preprint (2000) 

[28] Stroock, D.W. and Varadhan, S.R.S.: On the support of diffusion processes 
with applications to the strong maximum principle. In Proc. 6-th Berkeley 
Symp. Math. Stat. Prob., Vol III, Berkeley: Univ. California Press, 1972, pp. 
361-368 

[29] Wu, L.: Uniformly integrable operators and large deviations for Markov pro- 
cesses. J. Funct. Anal. 172, 301-376 (2000) 



20 



